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ABSTRACT 

The combination of large size, high stehar density, high metahicity, and Sersic surface brightness 
profile of the spheroidal component of the Andromeda galaxy (M31) within i^proj ~ 20 kpc suggest 
that it is unlike any subcomponent of the Milky Way. In this work we capitalize on our proximity to 
and external view of M31 to probe the kinematical properties of this "inner spheroid." We employ 
a Markov chain Monte Carlo (MCMC) analysis of resolved stellar kinematics from Keck^^/DEIMOS 
spectra of 5651 red giant branch stars to disentangle M31's inner spheroid from its stellar disk. We 
measure the mean velocity and dispersion of the spheroid in each of five spatial bins after accounting 
for a locally cold stellar disk as well as the Giant Southern Stream and associated tidal debris. For 
the first time, we detect significant spheroid rotation (v^ot ^ 50 km s~^) beyond i^proj ^ 5 kpc. The 
velocity dispersion decreases from about 140 km s~^ at R^to] = 7 kpc to 120 km s~^ at i^proj = 14 kpc, 
consistent to 2cr with existing measurements and models. We calculate the probability that a given 
star is a member of the spheroid and find that the spheroid has a significant presence throughout the 
spatial extent of our sample. Lastly, we show that the fiattening of the spheroid is due to velocity 
anisotropy in addition to rotation. Though this suggests that the inner spheroid of M31 more closely 
resembles an elliptical galaxy than a typical spiral galaxy bulge, it should be cautioned that our 
measurements are much farther out (2 — 14 reff) than for the comparison samples. 
Subject headings: galaxies: spiral — galaxies: kinematics and dynamics — galaxies: individual (M31) 
— techniques: spectroscopic — galaxies: local group 



1. INTRODUCTION 

The elegant progression of the Hubble sequence from 
ellipticals to spirals demonstrates that galaxy morphol- 
ogy can be described in large part based on the rel- 
ative importance of spheroid and disk subcomponents. 
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While it is now clear that the simple evolutionary path 
from elliptical ^^e arly-type " to disk-dominated "late- 
type" galaxies that Hubble ( 1936 ) originally proposed is 



incorrect, the physical origin of the Hubble sequence and 
the formation of and relationship between the different 
structural subcomponents remain subjects of vigorous re- 
search. 

The central spheroids of spiral galaxies fall into two 
categories, which can be expla ined by distinct forma- 
tion m echanisms as reviewed by |Kormendy fc Kennicutt 
(2004). Classical bulges, which are typically described 
as elliptical galaxy analogs with random stellar veloc- 
ity distributions, large velocity dispersions, and r^/^ de 
Vaucouleurs surface brightness profiles, are likely formed 
through violent merger/accretion events. Pseudobulges, 
which are more fiattened, have more ordered kinemat- 
ics, and have roughly exponential brightness cutoffs (or, 
more generally, Sersic profiles with low n-sersic 

values) 

are likely formed through secular heating of the disk. 
More detailed observations, yielding constraints on the 
structure and dynamics of bulges, will lead to a clearer 
understanding of possible formation scenarios. 

Any study of the inner regions of a galaxy is com- 
plicated by the presence of several spatially overlapping 
structural subcomponents, such as the disk, spheroid, 
and halo. Deconvolving these subcomponents to reveal 
the behavior of a single one is difficult. T raditionally, 
codes such as GALFIT ( |Peng et al.]|2002l ) or GIM2D 
(Simard |2002[ ) are employed to tit galactic integrated 
liffht profiles with the sum of a Sersic bulge and expo- 



nential disk (e.g. Courteaii||1996 Courteau et al.||201l[ ). 
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This technique is the only possible method for charac- 
terizing the structure of distant galaxies, but it suffers 
from strong assumptions about the characteristic light 
profiles of bulges and disks. In addition, degeneracy in 
the best-fit derived parameters can cloud interpretation 
of the results. 

Resolved stellar kinematics offer a complementary ap- 
proach to structural deconvolution of the nearest galax- 
ies. Instead of assuming specific surface brightness pro- 
files of disks and spheroids, one must only make the 
geometrical argument that a stable disk - a thin, flat 
structure - is kinematically colder (has a higher v^ot/cFy) 
than a stable spheroid. Separate components can then be 
identified and characterized by their distinct stellar ve- 
locity distributions. T he proximity of Andro meda (M31) 
at about 785 kpc (e.g.,|McConnachie et al.|2008) renders 
it the only large spiral galaxy other than the Milky Way 
(MW) where detailed photometric and kinematical ob- 
servations are possible with current observing facilities. 

We use resolved stellar kinematics to study M31's kine- 
matically hot "inner spheroid" at projected radii of 2-20 
kpc. Any description of this region - or that in the the 
intermediate region of any large galaxy - is necessar- 
ily complex; the literature is full of vocabulary such as 
"bulge," "spheroid," "inner spheroid," "outer spheroid," 
"disk," "thin disk," "thick disk," "extended disk," and so 
on. There is not yet a consensus on the best combination 
of these nouns to represent M31. For the purposes of this 
paper, we use the word "spheroid" to describe a kinemat- 
ically hot component: some combination of bulge, halo, 
and/or any other spheroidal component. Likewise, we re- 
fer to the kinematically colder population as the "disk," 
where this term includes any distinct disk components 
that may be present, such as the thin, thick or extended 
disks. 

Despite the possibility of multiple components, the 
disk is likely to be locally kinematically cold. 



Collins 



et al. (2011) claim that M31's stellar disk at rproj ^10- 
4U kpc consists of a cold thin disk and a warm thick disk, 
as is the case for the MW. Given that most of our fields 
are closer to M31's center than the innermost field of 



Collins et al. ( ^2011 ), and given their finding that the thin 
disk has twice the density of the thick disk and a shorter 
radial scale length, we expect the cold thin disk to dom- 
inate the stellar d isk population in our fields. Similarly, 



[Ibata et al.|( 2005| suggest that no more than about 10% 
of the totlil disk luminosity may be due to an extended 
disk component which also lags the cold disk. Though 
we do not know a priori the relativ e co ntributions of the 
thin, thick and extended disks, in §|4.2|we show that our 



assumption of a dominant thin component is justified for 
the purposes of measuring the kinematical parameters of 
the inner spheroid. 

Unlike the case of its stellar disk, M31's inner spheroid 
has no analog in the MW and is therefore of great in- 
terest. The spheroidal system at these ra dii in the MW 
irrelatively metal-poor {{Fe/H) ~ —1.6; iCarollo et al. 



[2007^ , is composed entirely of old stars, and has a power 
law spatial density profile of the form r^^p^^j , which cor- 
responds to an power-law surface brightness profile. 
Models such as t hose pro posed by Bullock & Johnston 
( |2005J and Zolotov et al.| ( 2010) suggest that the MW 
halo represents a population of accreted dwarf satellite 



galaxies. In contrast, the inner spheroid in M31 more 
closely resembles a bulge than a halo. It is more metal- 



enhanc ed than the MW halo, with [Fe/H] -0.7 (Kali- 



rai et al. 2006a), and has a Sersi c surface brightness pro- 
tile with nsersic ~ 2-4 (Pritchet fc van den Bergh||1994 
Guhathakurta et al. 2005 ; Courteau et al. 20lT| . In addi- 
tion, the stellar population of M31's spheroid is younger 
than that of the MW inner halo on average, with 4 0% of 

lia 



the stars younger than 10 Gyr (Brown et al. 



the stellar density is also significantly hi gher t 



2006k, and 



lan that at 



an equivalent location in the MW (Rei tzel et al.||1998) . 

The inner spheroid straddles territory between two 
well-studied components of the sphero id: the classical 
and bo xy bulges interior to ^ 1 kpc ([Athanassoula fc 
Beaton|[2QQ6| [Beaton et aL]|2007[ [Courteau et al.||2Ulip 



and the outer halo which dominates past 



(e.g. Guhathakurta et al. 



2005 



Irwin et al, 



2005 



5DTBc 



Ibata 



^et al. 2007) . In the central kpc, where the density is too 
high foFTeso lved stellar population spectroscopy, |Saglia| 



et al. (2010) analyzed integrated- light kinematics to re- 
veal a bulge rotation speed of 70 km s~^ and a velocity 
dispersion of 140 km s~^ at i^proj = 1.1 kpc on the major 
axis. However, they cautioned that this measurement is 
contaminated by the kinematically cold disk which may 
contribute nearly a third of the light at this radius. 

Farther out in M31's halo, kinematical surveys of 
the resolved stellar population using the Keck/DEIMOS 
multiobject spectrograph have mapped the cold sub- 
structure, as well as the underlying smooth virialized 
population, out to i^proi > 1^0 kpc (Guhathakurta e t al. 
2005; Chapman et al. 2006; Kalirai et al. 2006a; GilBerl 
et al. 2007n20Q9| ). [Chapman et al.| (2006) compiled kine- 
matics of ^ 1200 red giant branch (RGB) halo stars in 
scattered fields between R 



proj 



=8 and 70 kpc. Using 
a windowing technique to eliminate stars whose veloc- 
ities were consistent with that of the disk, they found 
that the velocity dispersion of the remaining popula- 
tion decreased radially outwards: cr^(i^proj) = (152 — 
0.9 i^proj kpc~^) km s~^. Subsequently, as part of 
the Spectroscopic and Panchromatic Lan dscape of An- 



drome da's Stellar Halo (SPLASH) survey, Gilbert et al. 
(2007) fit a double Gaussian profile to the velocity dis- 



tribution of RGB stars in a large contiguous region along 
the southeastern minor axis of the galaxy and measured 
a constant velocity dispersion of 128.9 km s~^ between 
^proj = 10 and 30 kpc. 

In recent years, the focus of SPLASH has migrated 
inwards, first to targe t the dwarf ga laxies Andromeda^^ 
I and Andromeda X 



_(|Tollerud'| 20"TT|) , N GC 205 (|Geha 
^al. 2008) ana^32 ([Howley et aT 



et al.| | 2006||Howley et _ ^ 

and now towards the disk- and bulge-dominatec 



inner regions of M31. The majority of the data ana- 
lyzed in the present paper come from the most crowded 
area targeted to date: a large contiguous disk-dominated 
area on the NE major axis with i^proj = 2-19 kpc. 
This area was selected to overlap the coverage of the 
Panchromatic Hubble Andromeda Treasury (PHAT) sur- 
vey, a five-year Hubble Space Telescope (HST) MultiCy- 
cle Treasury (MCT ) program that began in 2010 ( |Dal-| 
canton et aL]|2012[ ). 



I'he disk and spheroid share the inner regions of the 
galaxy with remnants of tidally disrupted galaxies. The 
dominant features in star-count maps of the 2-20 kpc 
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region are the Giant Southern Stream (GSS; the rem- 
nant of a tidahy stripped satehite galaxy) and the shelves 
(sharp edges in stellar den sity) created by it (Ibata et al. 
|2QQT| [Fardal et al][2QQ7) ). Both the GSS and a ^^sec- 



ondary stream," which is cospatial with the GSS but 
separated by 100 km s~^ in velocity, have been k inemat 



ically detected in multiple fields so uth of M31 (Kalirai 
et al.||2006b| [Gilbert et al.||2009D. [Fardal et al.] p)07| ) 



identified the northern extension of the stream in the 



Chapman et ah] ( |2006 ) and "Ibata et al." (2005) sample 
of RGBs and m the planetary nebulae of Merrett et al. 



( |2006D . 

'inhere are two principal challenges to a resolved stel- 
lar population kinematical study of the crowded inner 
spheroid of M31. First, we must select intended stellar 
targets whose spectra are least likely to be contaminated 
by close (in projection) stellar neighbors. Second, we 
must disentangle the stellar disk from the spheroid pop- 
ulation we wish to characterize. This is especially impor- 
tant in disk-dominated fields, where fewer than 20% of 
the stars may belong to the spheroid. Note that though 
a large part of this paper will focus on accounting for 
the disk contribution, our analysis method is designed to 
elucidate the nature of the inner spheroid rather than the 
disk. We do not attempt to make a statement here about 
the rotation curve of the stellar disk or the presence of 
a thick disk. We plan to analyze these components in a 
future paper. 

This paper is organized as follows. In §|2] we explain 
our target selection techniques, spectroscopic observa- 
tions and radial velocity extraction. In §[3] we describe 
our method for isolating and characterizing the spheroid 
velocity distribution in each of five spatial bins. In §[4] we 
discuss the implications of our results; finally, we sum- 
marize our findings in §[5j 

2. OBSERVATIONS 

Our data set for this project is a compilation of three 
sets of RGB spectra, two of which are presented here 
for the first time. A detailed technical description of the 
spectroscopic slitmask design and data reduction is given 
in Howley et al. (2012). In this section, we describe the 
target selection criteria for the different data sets and 
give an overview of the data acquisition and reduction 
metho ds c ommon to all the observations. 

In § |2.1[ we o utline the three data sets used in this 
paper. In § |2.2[ we describe the source catalo gs f rom 
which we select our spectroscopic targe ts. In § |2.3[ we 
explain our target selectio n cr iteria. In § |2.4[ we provide 
the observing details. In § |2.51 we give a ru ndown of the 
data reduction process. In § |2.6| and § |2.7[ we measure 
velocities of individual stars and determine the quality 



of those measurements, respectively. Finally, in § |2.8 
we discuss the detection and velocity measurement o: 
serendipitously detected stars. 

2.1. Data Sets 

The spatial coverage of our three data sets is shown 
in Figure [T] Each rectangular outline represents a 
Keck/DEIMOS slitmask, which covers approximately 
16' X 4' and yields 200-270 useful spectra. (The actual 
footprint of a DEIMOS slitmask is not perfectly rect- 
angular.) The SE data sets consists of three slitmasks 
oriented along the eastern minor axis of M31 (violet in 



Figure [T]). The M32 data set includes the five slitmasks 
covering the compact elliptical galaxy M32 directly south 
of M31 and one slitmask on the SW major axis (blue in 
Figure [T|. The SE and M32 data sets were observed 
during the 2007 and 2008 fall seasons. 

Our newest data set, from October 2010, covers a por- 
tion of M31's northeastern major axis spanning the pro- 
jected radial range 0.33°-1.38° or 2-19 kpc from the cen- 
ter of M31 in the plane of the disk (red and black in Fig- 
ure [T]). Five of these 15 slitmasks (black) were chosen 
to maximally overlap the regions of existing photometry 
from the first year of the PHAT program. 

The pre-imaging and reduction processes are identical 
for all three data sets. The primary difference in data 
acquisition is in the target selection: isolated sources are 
hand-selected for the SE and M32 slitmasks fro m a sin- 
gle mon ochromatic ground-based source catalog ( Howley] 
et al.|2 012), while we use an automated series of statisti- 
cal techniques as well as limited color information from 
the PHAT survey to choose targets for the MOT slit- 
masks. 

2.2. Source Catalogs 

All targets are chosen from an i'-band 2° x 2° 
CFHT/MegaCam mosaic centered on M31, obtained 
in November 2004. We run the software package 
DAOPHOT (StetsonI [T994| on the image to identify 



sources, fit PSFs, and produce a PSF-subtracted resid- 
ual image. The final catalog consists of nearly 2 million 
unique sources. 

For our 2010 Keck/DEIMOS observing run, we also 
had access to data from the first round of observations 
of the PHAT program. The data are organized into 
12f0 X 6f5 bricks; three half-bricks were available at the 
time of our slitmask design. From these data, we created 
lists of different stellar populations: metal-poor ([Fe/H] 
< -1.3), metal-intermediate (-1.3 < [Fe/H] < -0.7), 
and metal-rich ([Fe/H] > —0.7) RGB stars, and and hot, 
massive main sequence stars selected on the basis of SED 
fitting to six-filter HST photometry. Though we do not 
treat stars differently based on subcategory membership 
in the present paper, the kinematics of these stellar pop- 
ulations and their relationship to the kinematics of the 
RGB populations will be presented in a future work. 

2.3. Isolated Target Selection 

Not all the sources in our catalogs are equally good can- 
didates for multi-object spectroscopy; we prefer isolated 
targets, those whose spectra are least likely to be con- 
taminated by light from neighboring objects. We clas- 
sify this contamination as either crowding or blending. 
We define a crowded source as one that has at least one 
neighbor detected by DAOPHOT that is bright and close 
enough to potentially interfere with the spectrum of the 
source. In contrast, we define a blended source as one 
that is identified by DAOPHOT as a single source but 
for which visual inspection of the PSF-subtracted image 
indicates that more than one object may be present. 

Furthermore, we only target stars in the apparent mag- 
nitude range 2{) < i' < 22 for M32 and SE slitmasks 
(modified to 20 < < 21.5 for MOT slitmasks). We se- 
lect a bright limit of i' = 20 because the tip of the red gi- 
ant branch is at = 20.5 in M31 and so the MW contam- 
ination fraction increases significantly in brighter stars. 
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Fig. 1. Twenty-four Keck/DEIMOS multiobject slitmasks overlaid on the [Choi et al.| (|2002|) KPNO Burrell Schmidt B-band mosaic 
image of M31. Colored rectangles outline the slitmasks from the different data sets: (violet), M32 (blue), and MCT (red and black, 
respectively, for slitmasks for which target selection was based on only the ground-based CFHT/MegaCam photometry /astrometry catalog 
versus slitmasks for which target selection was based on a combination of the PHAT and CFHT/MegaCam photometry /astrometry 
catalogs). 



The faint-end limit is chosen because in very crowded ar- 
eas, it is difficult to recover high-quality spectra of stars 
fainter than about = 22. The surface density of sources 
in the area covered by the MCT slitmasks is so high that 
we can efficiently pack targets on our slitmasks even with 
a conservative faint-end limit of = 21.5. 

To choose targets for all slitmasks, we ffist sort possi- 
ble targets into three lists (1, 2, 3) in decreasing o rder o f 



isolation, with list 3 reserved for rejects (§ |2.3.1| - § |2.3.4 ). 
Within each list, we prioritize possible targets by magni- 
tude, giving highest priority to intermediate-brightness 
stars with 20.5 < i' < 21.0. In the final target selec- 
tion process , we e xhaust each list before moving on to 
the next. ^12.3.51 describes additional selection criteria 



specific to the five PHAT-based slitmasks. 



2.3.1. Crowding in the SE & M32 Data Sets 

In the M32 and SE data sets, we use a neighbor- 
rejection test to eliminate crowded sources. We reject 
any star with at least one neighbor with a sufficient com- 
bination of proximity and relative brightness, i.e., sat- 
isfies the following empirical criterion determined from 
visual inspection of the CFHT/MegaCam image: 



^nbr < ^tg 




3.0 



(1) 
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TABLE 1 

Keck/DEIMOS Multiobject Slitmask Exposures 



Mask 


Observation 


a [J2000] 


5 [J2000] 
(°/:-) 


P.A. 


^exp 


Seeing 


No. of 


No. of Usable 


No. of Usable 


Name 


Date (UT) 


(h:m:s) 


n 


(min) 


FWHM 


Slits 


Target Velocities 


Velocities of 








(Success Rate) 


Serendipitously Detected 
Stars 



M32_l^ 


2007 Nov 14 


00 


42 38.28 


+40 


51 34.0 


+160.0 


2x20 


O'.'S 


199 


188 


(94%) 


72 


M32_2 


2008 Aug 03 


00 


43 03.82 


+40 


55 07.7 


+70.0 


3x20 


0'.'6 


189 


166 


(88%) 


27 


iVl «J _«J 




00 


43 11.60 


+40 


52 34.7 


— 110.0 


3x20 


0^7 


203 


132 


(65%)^ 


10 


M32 4 




00 


42 13.87 


+40 


54 44.2 


+105.0 


3x20 


0^6 


165 


137 


(83%) 


119 


M32 5 


onns AiiP- 04 


00 


49 1 S 8S 


+40 




—75.0 


3x20 


0^6 


177 


152 


(86%) 


72 




onns AiiP- S1 


00 


41 20.41 


-k41 
-|-^ J. 


51 32.2 


0.0 


■^y on _L 1x10 


0^7 


169 


152 


(90%) 


128 


SET 


2008 Sept 01 


00 


43 38.74 


+41 


10 17.4 


+39.0 


2x10 


0(^8 


170 


148 


(87%) 


50 


SE8 


2008 Sept 30 


00 


44 00.82 


+41 


09 27.1 


-113.0 


3x15 


0(^5 


197 


178 


(90%) 


27 




zUUo kJcXj Ui 


uu 


A A /in 0<5 


1 A 1 

+41 




— bU.U 


Z X iz.o + i X io 


n" A 
U. 4 


zU4 


ioO 


(91 /oj 


1 1 
ii 


mctA5 


2010 Oct 07 


00 


44 18.33 


+41 


39 28.1 


+270.0 


3x16 


0(^6 


212 


197 


(93%) 


82 


mctB4 


2010 Oct 08 


00 


44 29.04 


+41 


35 10.0 


+270.0 


3x16 


0(^7 


177 


172 


(97%) 


81 


mctC3 


2010 Oct 07 


00 


45 11.89 


+41 


53 37.7 


+90.0 


3x16 


0(^5 


198 


170 


(86%) 


69 


mctD3 


2010 Oct 07 


00 


45 09.46 


+41 


49 08.8 


+90.0 


3x18 


0^/6 


209 


185 


(84%) 


69 


mctE3 


2010 Oct 08 


00 


46 53.23 


+42 


14 59.3 


+90.0 


3x17 


0^/7 


221 


202 


(91%) 


18 


mct04p 


2010 Oct 07 


00 


44 51.81 


+41 


25 19.2 


-142.3 


3x16 


0(^6 


254 


223 


(88%) 


24 


mctOSp 


2010 Oct 08 


00 


44 19.70 


+41 


32 53.5 


-52.3 


3x16 


0(^5 


254 


188 


(74%) 


100 


mctOGp 


2010 Oct 07 


00 


44 33.71 


+41 


36 17.6 


-52.3 


3x16 


0(^5 


251 


211 


(84%) 


43 


mctOTp 


2010 Oct 08 


00 


44 41.77 


+41 


40 03.0 


-52.3 


3x16 


0(^6 


264 


210 


(80%) 


71 


mct09p 


2010 Oct 08 


00 


45 39.23 


+41 


38 39.1 


-142.3 


3x16 


0(^6 


252 


213 


(85%) 


22 


met 1 Op 


2010 Oct 07 


00 


45 08.24 


+41 


46 19.2 


-52.3 


3x18 


0(^9 


255 


207 


(82%) 


34 


mctl2p 


2010 Oct 07 


00 


45 28.34 


+41 


53 23.3 


-52.3 


3x18 


0(^9 


265 


212 


(80%) 


12 


mctl3p 


2010 Oct 08 


00 


45 42.02 


+41 


56 42.4 


-52.3 


3x18 


0(^7 


259 


217 


(84%) 


23 


mctlSp 


2010 Oct 08 


00 


45 54.36 


+41 


59 43.1 


-52.3 


3x17 


0(^9 


261 


206 


(79%) 


10 


mctlGp 


2010 Oct 08 


00 


46 08.44 


+41 


02 58.6 


-52.3 


3x18 


0(^7 


258 


221 


(86%) 


5 



Total: 



5263 



4472 (85%) 



1179 



^ The "M32_l" slitmask was originally named "M32" at the time of submission of the slitmask design. 

^ Due to a warp in the slitmask during the time of observation, approximately 30% of the slitlets from mask M32_3 did not produce useful spectra. 



Here, /tgt and /nbr are the i'-band magnitude of the tar- 
get source and the neighbor, respectively, and d is the 
distance in arcseconds between the target and the neigh- 
bor. This cut ehminat es about 90% of th e stars in the 
M32 and SE data sets ( |Howley et aL]|2Q12 ). 



2.3.2. Blending in the SE & M32 Data Sets 

We identify hkely blends in the SE and M32 data sets 
by visually inspecting the high-pass filtered and PSF- 
subtracted versions of the z'-band CFHT/MegaCam im- 
age at the locations of the stars that survive the crowd- 
ing test. Each target is flagged as unblended, marginally 
blended, or badly blended depending on the degree to 
which its image resembles the PSF on the high-pass fil- 
tered image and the strength of systematic residuals at 
its location on the residual image. 

2.3.3. Crowding in the MCT Data Set 

Based on our experience with the M32 and SE spec- 
troscopic data sets, we decide to use a slightly relaxed 
crowding criterion for the MCT data set, rejecting (as- 
signing to list 3) any catalog entry with at least one 
neighbor which satisfies the following: 



^nbr < ^1 



tgt 




(2) 



This change accounts for the fact that the seeing during 
our spectroscopic observations tends to be slightly better 



than the 0^8 CFHT seeing upon which the original cri- 
terion was based. This cut eliminates 80% — 90% of the 
stars in the inner slitmasks, 50% — 70% at intermediate 
radii, and only 40% in the furthest slitmask, mctE3. 

2.3.4. Blending in the MCT Data Set 

In order to avoid a tedious visual inspection of the 
large area covered by the MCT data set, we design 
two empirically-based statistical tests to detect possible 
blends. We visually inspect and flag as "blended" or 
"non-blended" < 100 objects in each of three small rep- 
resentative image sections at different distances from the 
center of M31. We then design quantitative tests that 
approximately reproduce our visual classifications. 

The first test is based on the DAOPHOT-generated 
goodness-of-fit parameter chi and shape parameter 
sharp. Objects that appear isolated based on visual in- 
spection fall into a well-defined locus in chi/sharp space, 
as shown in Figure [2] Based on this relationship, we re- 
tain only objects with sharp < 0.2. We assign a radially 
dependent linear cut in chi, accepting all stars with chi 
< 0.3 in the crowded areas and progressing to the more 
stringent criterion chi < 0.2 in the least crowded outer- 
most slitmask. This cut eliminates 20% of the remaining 
candidate targets. Because the dividing line between iso- 
lated sources and likely blends is less well defined in chi 
than sharp, we allow a buffer zone 0.5 units wide in chi. 
The 1% of stars in the buffer zone are relegated to list 2, 
but not rejected outright. 

In the second test for possible blends, we compare the 
apparent quality of subtraction to the normalized RMS 
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flux of a 5 X 5 pixel square of the PSF-subtracted im- 
age centered on the source. This value tends to increase 
with apparent degree of blending. We determine that the 
best cut is a linear function of magnitude, where blended 
sources have > 0.3 at = 20 and > 1.2 at 

<nux> <nux> 

= 21.5 (Figure [3|. The 36% of stars that fail this test 
are flagged as "possibly isolated" and pushed to list 2. 
Because the correlation between apparent PSF subtrac- 
tion quality and RMS is not as tight as those in the 
chi/sharp test, we do not use the RMS cut to reject 
(assign to list 3) targets that are "isolated" according to 
both the neighbor-rejection and chi/sharp tests. 

2.3.5. Target Selection for PHAT-based Masks 

We also design five slitmasks based jointly on CFHT 
data and the PHAT survey-based stellar population lists 
described in 



2.2 



We require that a star from the PHAT 
catalog be in the CFHT catalog and have passed through 
the filters described above to be considered for selec- 
tion. To ensure final selection of the most isolated (list 1) 
PHAT objects, we push all non-PHAT objects down one 
list. For these five slitmasks, then, the priority scheme is 
as follows: list 1 consists of only isolated, PHAT-selected 
sources; list 2 includes the isolated CFHT-only sources; 
and list 3 includes all of the possibly-isolated sources. 
Because the shape of the PHAT survey bricks is differ- 
ent from the shape of DEIMOS slitmasks, large areas of 
the DEIMOS slitmasks that overlap PHAT survey bricks 
do not target PHAT-based sources. In these areas, the 
mask design software automatically proceeds to list 2 to 
select CFHT-based objects. 

2.3.6. Summary of Target Selection 

To summarize, each source passes through three tests 
for isolation: neighbor-rejection, chi/sharp, and RMS. 
Each source receives a score for each test: if isolated, 
0.6 if marginally isolated, and 2.0 if not isolated. The 
sum of the three scores determines which list the source 
belongs in. A total score of maps to list 1; a total of 
less than 2.0 maps to list 2; and a star with a total score 
of 2.0 or greater is assigned to list 3. Hence, any star 
that fails either the neighbor-rejection or chi/sharp test 
cannot be selected, and any star that passes all three 
filters is given highest priority. 

None of the methods described here can fully eliminate 
the possiblity of placing one slit over several objects. Es- 
pecially in the crowded inner areas, it is very common to 
obtain mutiple spectra in one slit. We discuss our han- 
dling of these serendipitous detections in § |2.8[ In addi- 
tion, a small percentage of the target spectra may still 
be contaminated with light from nearby objects; objects 
with unusable spectra are identified by eye and removed 
from the sample at t he e nd of the data reduction pro- 
cess as described in 



2.7 



The target selection process 
outlined above serves simply to make educated guesses 
about the objects best suited for spectroscopy. 

In § |4.6| we show that the spectroscopic target selection 
criteria (PHAT CMD vs. magnitude only) and actual de- 
gree of crowding (isolated vs. sharing a slit with another 
bright object) have minimal, if any, effect on the mea- 
sured velocity distribution. 
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Fig. 2. — chi/sharp isolation criteria. Visually identified blends 
(red crosses) and non-blends (green circles) are shown from a small 
representative area close to the center of M31. We reject all sources 
with DAOPHOT goodness-of-fit parameters chi > 0.3 (to the right 
of the dashed line) or shape parameter sharp > 0.2 (above the 
dotted line). The chi cutoff value is lowered in less dense fields 
farther from the galactic center. 
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Fig. 3. — RMS isolation criterion. Visually identified blends (red 
crosses) and non-blends (green circles) are shown from a small rep- 
resentative area about 8 kpc from the center of M31. The RMS 
deviation of the pixel values of the residual image in a 5 x 5 pixel 
square centered on each source is plotted against the i' magnitude 
of the source. We reject all sources above the solid black line. 

2.4. Observations 

All slitmasks are observed using Keck/DEIMOS with 
the 1200 line mm~^ grating. This configuration yields a 
spatial scale of 0^12 pixel" ^ and a spectral dispersion of 
0.33 A pixel"^. We set the central wavelength to 7800 
A, corresponding to a wavelength range of ~ 6450 — 9150 
A. The exact wavelength range for each slit varies as a 
result of location on the slitmask and/or truncation due 
to vignetting. The wavelength region is chosen to tar- 
get the Ca II triplet absorption feature present in RGB 
stars. The anamorphic distortion factor for this grat- 
ing and central wavelength is 0.606. Therefore, each 0V8 
wide slitlet subtends 4.1 pixels. Better still, excellent see- 
ing conditions (~ 0^6) during observations can provide 
somewhat better spectral resolution yielding an average 
resolution of 3.1 pixels = 1.0 A. 

Reliable spectra (those that yield secure velocities, as 
described in §|2.7|), are obtained from 4465 of the 5263 
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slitlets. Approximately 30% of the shtlets from shtmask 
M32_3 did not produce useful spectra due to a warp in 
the slitmask. 

See Table [l] for information on the positions of the slit- 
masks and the number of useful spectra recovered from 
each one. 

2.5. Data Reduction 

The Keck/DEIMOS multiobject slitmasks are pro- 
cessed using the spec2d and specld software (version 
1.1.4) developed by the DEEP Galaxy Redshift Survey 
team at the University of California, Berkeley (Davis 
et al. 2003 )q| Briefly, the reduction pipeline rectifies, 
tlat-tield anairinge corrects, wavelength calibrates, sky 
subtracts, and cosmic ray cleans the two-dimensional 
spectra, and ext racts the one-dimens ional spectra. For 
more details, see Howley et al. (2012). 



2.6. Cross- Correlation Analysis 

Line-of-sight (LOS) velocities for resolved target s are 
measu r ed fro m the one-dimensional spectra using a Geha 
et al. (2010) modified version of the visual inspection 
software zspec^ developed by D. Madgwick for the DEEP 
Galaxy Redshift Survey at the University of Califor- 
nia, Berkeley. The software determines the best-fit 
LOS velocity for a target by cross-correlating its one- 
dimensional science spectrum with high signal-to-noise 
stellar templates in pixel space and locating the best fit 
in reduced-x^ space. 

A-band telluric corrections and heliocentric corrections 
are calculated and applied to the measured LOS veloci- 
ties. The A-band telluric corrections, which account for 
velocity errors associated with the slight mis-centering 
of a star in a slit, are deter mine d usin^ the method di s- 
cussed in |Sohn et ah] ( |2007 l) and iSimon fc Geha| ( [20071 ) . 

LOS velocity errors are determined for each star by 
adding in quadrature the cross-correlation based velocity 
error and a systematic error estimated by repeat obser- 
vations, as described in Howley et al. ( 2012| ). The typical 
LOS velocity error in our sample is 4-5 km s~^. 

2.7. Quality Assessment 

Each two-dimensional spectrum, one-dimensional 
spectrum, and corresponding Doppler shifted template 
match are visually inspected in zspec and assigned a qual- 
ity code based on the reliability of the fit. This process 
allows the user to evaluate the quality of a spectrum 
and reject instrumental failures and poor quality spec- 
tra. Velocity measurements based on two or more strong 
spectral features are labeled "secure." Velocity mea- 
surements based on one strong feature plus additional 
marginal features are labeled "marginal." Spectra that 
contain no strong features, low S/N and/or instrumen- 
tal failures are considered unreliable, and so are not in- 
cluded in our analysis. Addi tional details on quali ty code 
assignment can be found in Guhathakurta et al. (12006). 
During this process, we also identify and Hag 43 likely 
MW M dwarfs based on their strong surface-gravity sen- 
sitive Na I 8 190 A doublet. These are excluded from the 
radial velocity analysis. 



2.8. Serendipitous Sources 

Upon visual inspection of the one-dimensional and two- 
dimensional spectra during the quality assessment phase 
outlined in § |2.7[ some fraction of the slits clearly show 
that the slitlet intersects more than one star: the target 
star and one or more serendipitously detected stars, or 
serendips. Serendips are detected via one of two meth- 
ods: through continuum detections that are offset from 
the primary target in the spatial direction, or by the 
detection of spectral features that are offset from the 
primary target in the spectral direction. Serendip de- 
tections occur frequently in the inner parts of M31 and 
close to M32 due to the severe crowding and blending in 
the CFHT/MegaCam data. Serendips are also assigned 
quality codes and those with secure or marginal veloc- 
ities are included in the radial velocity analysis. More 
information on seren dipitous detections can be found in 
Howley et aL](|2012[). 



3. DATA ANALYSIS 

We perform our analysis in each of the five spatial re- 
gions labeled in Figure |4j These regions are the SE minor 
axis (not expected to yield constraints on the spheroid 
rotation velocity); the SSW quadrant (expected to yield 
a constraint on the rotation velocity via a negative veloc- 
ity offset from the systemic velocity of M31); and three 
regions along the northeast major axis, which we name 
NEl, NE2 and NE3 in order of increasing projected ra- 
dial distance (expected to yield three independent esti- 
mates of the rotation velocity via positive velocity offsets 
from the systemic velocity of M31). Note that these re- 
gions are defined by lines of constant position angle and 
projected radius, and are not quite the same as the three 
data sets shown in Figure [l] The positions of and num- 
ber of stars within each region are given in the first six 
columns of Table |2j 

The kinematically cold peak at around —100 km s~^ in 
Figure |5] suggests that our measured velocity distribution 
has a significant contribution from the disk. Hence, it is 
imperative that we realistically account for disk contribu- 
tions. Instead of adopting a specific model for its velocity 
field, we only assume that the stellar disk is locally cold 
with a symmetric velocity distribution. As explained in 



3.1 , we apply this assumption to divide each region into 



several subregions. To each subregion we fit two Gaus- 
sian distributions, corresponding to a kinematically cold 
and a hot component, where the hot component is re- 
quired to have the same mean velocity and velocity dis- 
persion across all subregions in a region (described in 



3.2). Lastly, in §3.3 we describe how we modify our 



http: //astro. berkeley.edu/""cooper/deep/spec2d/ 



analysis to account for the possibility of contamination 
by the CSS and associated tidal debris. 

3.1. Choice of Subregions and Expected Disk LOS 
Velocity Pattern 

Because we assume that the disk is only locally cold, we 
fit for the disk in each of many small subregions. The spa- 
tial boundaries of individual subregions are dictated by 
two competing desirable factors, namely a small spread 
in disk mean velocity and high number statistics. Our 
assumption of a perfectly cold disk is only strictly true 
in the limit of infinitely small subregions. On the other 
hand, a multi-Gaussian fit to a velocity distribution re- 
quires a somewhat large number of points. We arbitrarily 
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decide that 100 is the minimum satisfactory number of 
points. 

To estimate the spread of the mean disk velocity as 
a function of position, we employ a simple geometrical 
model for the rotation pattern of an inclined disk with 
perfectly circular motion ( (Guhathakurta et aL,,1988| ): 



Proiected distance from center of M3 1 (kpc) 
10 5 



^rot(^deproj)sin(z) 

V'l + tan2(APA)/cos2(z) 



(3) 



where ^, r] are tangent-plane coordinates with origin at 
the center of M31, i = 77° is the inclination of the disk 
of M31, Vsys = —300 km s~^ is the systemic heliocen- 
tric velocity of M31, APA is the position angle projected 
onto the plane of the sky measured relative to the ma- 
jor axis of the disk of M31, rdeproj is the radial position 
measured in the plane of the disk, Vj-ot is the disk ro- 
tation speed, and the -|- and — signs apply to the NE 
and SW halves of the disk, respectively. An azimuthally 
averaged estimate for '^rotl^deproj), based on HI kinemat- 
ics, ranges from ~ 250 km at rdeproj = 15 kpc to 
- 175 km s-i at 5 kpc (ICorbem et al.||2010[ ). The ex- 
pected velocity spread calculated in this way will not be 
exact, because in addition to APA, the true spread in 
mean stellar velocity over a subregion is influenced by 
Ardeproj of the subregion, any departure from perfectly 
circular motion, and the intrinsic local velocity distribu- 
tion (due to a multiple-component disk, for example). 

Using Equation [3] with Vrot = 250 km s~^, we bin 
our data in each region into subregions based on po- 
sition angle. The angle subtended by a single subre- 
gion is approximately the greater of two APA crite- 
ria: 1) the APA such that the change in Vobs over a 
subregion due to APA is 10 km s~^, or 2) the APA 
that includes 100 data points. Our final subregions 
are shown in Figure Hj In the NE1-NE3 and SSW 
regions, we identify these subregions with subscripts 
that increase with distance from the nearest major axis: 
NEli,NEl2,...,NEl7; SSWi, . . . , SSW5, etc. In the SE 
region, the outer, inner south, and inner north subre- 
gions are named SEi,SE2,SE3, respectively. Note that 
we use this rotation pattern only to estimate appropri- 
ate bin sizes, not to determine disk rotation speeds. The 
final positions of and number of stars in each subregion 
are presented in the first six columns of Table |4] in the 
Appendix. 

3.2. Fitting the Velocity Distribution Model 

We use a Markov chain Monte Carlo (MCMC) sam- 
pler to find the velocity distributions of the disk and 
spheroid in each region. The spheroid is modeled as a 
single Gaussian distribution for each region, while the 
disk is modeled as a Gaussian distribution for each sub- 
region. The model likelihood at each point in parameter 
space is then: 



N 



'j I ^sph,r(j); '^sph,r(j) J 



(1 - fr(j)) J^{vj I Vdisk,s(i)> <^disk,s(j)) (4) 




1.0 0.8 0.6 0.4 0.2 
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Fig. 4. — Division into regions and subregions. Each point cor- 
responds to a single velocity measurement. Solid lines delineate 
the five regions for which we obtain independent kinematical pa- 
rameters of the spheroid: inner northeastern (NEl), intermediate 
northeastern (NE2), outer northeastern (NE3), southeastern (SE), 
and south-southwestern (SSW). Colors show the subregions within 
each region that are used to determine the disk contribution to 
the velocity distribution. With the exception of the SE region, our 
subregion naming convention is based on distance from the major 
axis. For example, in the NE3 region, NE3i (light gray) strad- 
dles the major axis, NE32 (black) is slightly farther out, and NE33 
(dark gray) is farthest from the major axis. In the SE region, we 
use projected distance from the center of M31 in addition to lines 
of constant position angle to define three subregions: SEi (outer, 
yellow), SEI2 (inner south, dark brown) and SEI3 (inner north, 
light brown). 



300 




where N is the number of stars in our data set, and 
star j with measured velocity Vj is found in region r(j) 



-400 -200 

Stellar heliocentric velocity (km/s) 

Fig. 5. — Distribution of heliocentric radial velocities for all 3697 
stars with reliable velocity measurements from the MCT slitmasks. 
The peak around —100 km s"-*^ corresponds to the disk, the LOS 
projection of whose rotation speed causes stars to be offset by as 
much as +200 km s"-*^ with respect to the —300 km s"-*^ systemic 
velocity of M31 (dotted hue). 

and subregion s{j). The notation J\f{x | /i, a) indicates 
the Gaussian distribution with mean ji and standard 
deviation <j evaluated at x. The scalar /g, subject to 
< /s < 1, is the fraction of the stars in subregion s 
that belong to the spheroid. The spheroid in region r 
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has mean velocity Vsph,r and velocity dispersion crsph,r, 
while the disk in subregion s is characterized by velocity 
^disk,s and dispersion crdisk,^. 

Subregion SSW4 includes stars from the galaxy M32, 
so for that subregion only we replace the single- Gaussian 
disk model A/'(i'j | ^disk,s(j)7 <^disk,s(j)) with a double- 
Gaussian model. 

The likelihood in Equation [4] can be sampled indepen- 
dently in each of the fi ve regions. We use the code emcee 
(|Foreman-Mackey et al.||2Q12]) , wh ich implements the 
afiine-invariant ensemble sampler of [Go o dman &: Weare| 
([2010 ), to perform the MCMC algorithm. In addition 
to the likelihood above, we must specify prior probabil- 
ity distributions for the parameters. For fs we use a 
uniform distribution on the interval [0, 1]; for the mean 
velocities ^sph,r and ^disk,s we use flat priors; and for the 
dispersions (Jsph,r and crdisk,s we demand positive val- 
ues. If we sample each region independently, then the 
parameter space for region r includes ^sph,r, crsph,r, plus 
the set of /g, ^disk,^, crdisk,s for each subregion s within 
r. We initialize the MCMC at ^sph,r = -300 km s"\ 

c^sph,r = 150 km s"\ each fs = \, and ^disk,s and crdisk,s 
= the sample mean and standard deviation of the veloc- 
ities in subregion s. 

Briefly, the MCMC ensemble sampler works as follows. 
It explores the parameter space by maintaining a set of 
walkers. Each walker represents a point in the parame- 
ter space. At each iteration of the MCMC, each walker 
takes a step in the parameter space by choosing another 
walker and stepping along the line in parameter space 
connecting itself to the other walker. The step size is 
chosen stochastically and allows interpolation as well as 
extrapolaton. In effect, the walkers choose their steps 
based on the covariance of the set of walkers. After each 
step is taken, the posterior probability distribution at 
the new point in parameter space is evaluated. Steps 
that increase the probability are always accepted, while 
steps that decrease the probability are sometimes ac- 
cepted. After a large number of steps, the ensemble of 
walkers will sample the parameter space with frequency 
proportional to the posterior probability distribution; we 
can draw fair samples from the distribution by selecting 
points from the histories of the walkers. We estimate 
the mean and variance of each parameter (assuming uni- 
modal distributions) based on sample statistics of the 
histories of the walkers. In particular, we allow each of 
32 walkers to take 10,000 steps. We then compute the 
mean spheroid velocity as the mean value of ^sph,r over 
the last 2000 steps of all the walkers (i.e., the mean of 
64,000 points), and the 68% confidence interval as the 
standard deviation of that quantity over the last 2000 
steps of all of the walkers. We use the same method 
to estimate the values and uncertainties of the spheroid 
dispersion and the disk parameters. We report these val- 
ues in columns 7 and 8 of Table [2] and illustrate them in 
Figure |9] 

3.3. Accounting for Tidal Debris Associated with the 
Giant Southern Stream 

The analysis described thus far separates the spheroid 
from the disk (and from M32 in subregion SSW4). The 
mean velocity and dispersion we extract describe the av- 
erage properties of the spheroid, regardless of its under- 



lying structure. 

Several lines of evidence, however, suggest that our 
measurements do not well represent a rotation curve of a 
smooth spheroid. First, the mean spheroid velocities in 
the five regions do not follow a physical rotation pattern: 
the mean velocity of the NE2 region is consistent with 
zero, while the surrounding regions (NEl and NE3) are 
rotating at about 50 km s~^ in the same direction as the 
disk. Second, the sum of hot and cold Gaussians does not 
well represent the velocity distributions. A chi-squared 
analysis applied to the data binned by 20 km s~^ re- 
veals that that the probability of the model representing 
the velocity distribution is relatively low (see the final 
column of Table [2| . 

Velocity histograms, such as those in Figure [6j sug- 
gest that cold substructure, possibly tidal debris from 
the GSS, could be skewing our measurements. A close-up 
view of the negative-velocity tail of the NE2 histogram 
reveals a cold spike of about 10 stars in excess of the 
Gaussian tail at —580 km s~^ (Figure |6j top right). A 
slight overdensity of points at this velocity can be seen 
in the NEl and NE3 histograms as well. Though noth- 
ing is immediately visible in the SE or SSW regions, a 
few extra stars around —580 km s~^ would be partially 
concealed by the bulk of the velocity distribution. 

Figure [7] shows another projection of these data: the 
velocity 01 each star in the NEl, NE2, and NE3 regions 
plotted against i^proj- Also shown (in turquoise trian- 
^les) are data from , two GSS fields at 17 and 21 kpc 
( Gilbert et al.||200"9 ). The concentrations of stars at ap- 
proximately —500 and —400 km s~^ in these data rep- 
resent the GSS and the secondary stream, respectively. 
The magnitudes of the central stream velocities increase 
with decreasing i^proj as the streams fall into the poten- 
tial well of M31 from the south. The black crosses in 
Figure [7| show the expected stream v elocity as a function 
of radius closer to the center of M31 ( Fardal et al.|2012 ). 
The NE region data from this work show clear concen- 
tr at ions of objects near the predicted stream velocity, 
continuing the trend seen south of the galaxy. Hence, it 
seems unlikely that the peak in the Figure [6] histograms 
is simply a binning artifact, and probable that it comes 
from the northern extension of a cold tidal stream. 

To account for the presence of the GSS and its asso- 
ciated tidal debris, we repeat the MCMC fits in regions 
NEl, NE2, NE3, and SSW after removing ah stars within 
cr = ±30 km s~^ of the predicted velocities of the two 
strea ms. (The measured v elocity dispersion of the stream 
from Gilbert et al. (20091) is 20 km s~^; however, we use 
the larger value to ensure that we exclude all stream 
stars despite the slightly uncertain mean stream veloc- 
ity.) In the SE minor axis region, we exclude all stars 
with velocity v < —600 km s~^ or > km s~^. The 
former is to account for GSS debris, and the latter for 
SE Shelf stars (GSS debris from third pericentric pas- 
sage) (Gilbert et al.| |2007). 

We account lor the fact that we have removed stars 
within a range of velocities by renormalizing the Gaus- 
sian distributions. That is, given a Gaussian velocity 
distribution, we compute the fraction of the probability 
mass that falls within the excised velocity range and scale 
up the remainder of the distribution so that the integral 
over the remaining (unexcised) velocities is unity. The 
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resulting kinematical spheroid parameters are reported 
in columns 7 and 8 of Table |3l 

The resulting fit for region NEl is plotted in Figure [Sj 
The first seven panels show the seven subregions in NEl. 
In each panel, each walker at the end of 10,000 steps is 
represented by a set of colored lines: a red Gaussian with 
parameters /g, ^sph,r, crsph,r and a blue or green Gaus- 
sian with parameters (1 — /g), ^disk,s7 crdisk,s- The sum 
of these two is shown in violet. These distributions are 
overlaid on the velocity histogram of stars in the subre- 
gion. The velocity ranges excised for stream debris are 
shown in two ways: the velocity range is shaded in light 
gray, and the stars in this range are colored red. 

The bottom middle panel shows the cumulative best-fit 
distributions: the red line is the Gaussian corresponding 
to the mean (^sph,r, crsph,r)- The blue curve is the cumu- 
lative disk distribution: the sum of the mean disk Gaus- 
sians from the seven subregions. The violet is the sum 
of the other two curves. The excised velocity ranges are 
respresented in the same way as in the subregion panels. 

The bottom right panel shows the positions of the 
walkers in parameter space after 10,000 steps. There 
are eight concentrations of points corresponding to the 
eight (^, ay) pairs (one spheroid and seven disk). Each 
pair is summarized by an ellipse displaying the mean and 
dispersion of that distribution. 

4. DISCUSSION 

This section is organized as follows. We first present 
the rotation curve and veloci ty d ispersion profile for the 



inner spheroid in §4.1 In § 4.2 we compare the veloc 



ity dispersion of the cold component to previous mea- 
surements of the stell ar di sk as a sanity check on our 
analysis method. In §|4.3| we show that the exclusion 



of velocity ranges corresponding to tidal debris from the 
GSS significantly impacts the measured kinematical pa- 
rameters of the spher oid, but has minimal effect on the 
cold component. In §|4.4| we compute a spheroid mem- 



bership probability for each star in our sample. In § 4.5 



we explain that the spheroid is likely supported by ve- 



locity anisotropy in addition to rotation. Finally, in § 4.6 
we show that neither spectroscopic target selection cri- 
teria nor degree of crowding introduces a significant bias 
towards the kinematically cold or hot population. 

4.1. Kinematical Parameters of the Inner Spheroid 

The spheroid distributions corresponding to the best- 
fit ^sph and (Jsph (hereafter Cy) are plotted in Figure [6] 
for each of the five regions. (This figure is simply a com- 
pilation of the summary panels of Figures Hand 13-16.) 
Kinematical and goodness-of-fit parameters, accounting 
for the GSS and associated tidal debris, are reported 
in Table [s] The probabilities have improved signif- 
icantly from the uncorrected values in Table [2] Four of 
the five regions now have probabilities greater than 90%, 
and that of the fifth (SSW) has increased by a factor of 
15. 

The dispersion profile is shown in the right panel of 
Figure [To) The profile appears to decrease smoothly with 
radius, though it is consistent with flat to 2a. The best- 
fit line to this profile is 

cr^(aeff) = (159.5±10.8)-(2.9±1.3)^^ km s-^ (5) 



where aeff is the effective spheroi d major axis coordi- 
nate, assum ing a 5:3 axis ratio (Pritchet & van den 
Bergh 1994). Our dispersion profile is consistent with 
that measured by Gilbert et al. (2007) but is slightly 
offset from other e xistmg measurements, including the 
Saglia et al. (2010) integrated-light measurement at 1.1 
kpc on the major axis of the bulge (blac k square in Fig- 
urep^ and the linear dispersion profile of [Chapman et al. 
( 20(IB") . These differences should be taken lightly, though, 
because the measurements are not directly comparable. 
The Saglia et al. (2010) point may be slightly defiated by 
contributions from the low-dispersion disk stars, which 
those authors est imate contribute a bout 30% of the light 
in the slit. The Chapman et al. (2006) profile, mean- 
while, is the innermost limit of measurements primarily 
made farther out in the halo, so we do not necessarily 
expect agreement at the radii covered by our study. 

We also compare our results to the dispersion pro- 
file produced by a model with spherical, isotropic, non- 
rotating bulge and halo stellar components as well as a 
stellar disk and halo {M/L = 2. 5) and bulge (M/L = 
5.6) (dotted orange line in Figure [iQl lFardal et al.|2012l ). 
The model falls below our data at tHe 2cr level in the iNE2 
and NE3 regions and at the 3a level in the SE and NEl 
regions. The mismatch suggests that, for example, the 
mass profile of the model galaxy is too shallow, or the 
gradient of the density profile of the tracer population is 
too large. A more detailed analysis is beyond the scope 
of this paper. 

Extracting the intrinsic rotation curve of the spheroid 
is nontrivial. The relationship between the mean LOS ve- 
locity and the intrinsic rotation velocity depends on the 
orbital dynamics of the spheroid, and in general is diffi- 
cult to determine without detailed 2D or 3D kinematical 
mapping. Therefore, in the left panel of Figure [To] we 
simply present the mean LOS component of the veloc- 
ity versus the effective major axis coordinate (based on 
a 5:3 axis ratio). Though two points are consistent with 
zero to 2(7, we detect significant rotation in the SSW, 
NEl and NE3 regions. The average value of \v — vm3i|, 
excluding the SE minor axis point, is 52.6 ± 6.8 km s~^ 
(solid line in Figure 10). All four off-minor axis points 
are consistent with this mean velocity to better than la. 
This is the first measurement of significant rotation in 
the inner spheroid. 

The velocity dispersion of M31's inner spheroid is sim- 
ilar to that the halo of the MW at similar radii, but 
its mean velocity is significantly larger. The velocity 
ellipsoid of the MW's inner halo is (cry^^ , cry^ , cry^ ) = 
(150 ± 2, 95 ±2, 85 ±1) (ICarollo et al.|2010|), on the same 
order as the LOS component of the M31 spheroid disper- 
sion. However, the MW's inn er halo has a mean r otation 
velocity consistent with zero dCarollo et aL][2QTQl ). 



4.2. Dispersion and Velocity of the Cold Population 

The cold component (^, ay) for each subregion are plot- 
ted against rdeproj in Figure 17 As a sanity check on our 
analysis method of fitting multiple Gaussians to subre- 
gions, we compare these values to previously measured 
kinematical parameters of the stellar disk. The cold 
component in a single subregion has an average disper- 
sion o f 58 km s~^, reasonably consistent with the Collins 
et al. ( [2011} dispersion measurements of 51 km s~^ and 
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Fig. 6. — Maximum likelihood fits of a kinematically hot spheroid to each of the five regions in our sample, after excluding the velocity 
range encompassing the GSS and its associated tidal debris (shaded pink). Violet lines show the cumulative region cold component; blue 
show the best-fit spheroid Gaussian; green show the sum of these two components. The dashed lines show the systemic velocity of M31 
relative to the MW. Individual subregion cold components are shown in orange in the SSW region panel, but left out of the other panels 
for clarity. We remind the reader that our fitting procedure makes use of the distinct velocity profiles of the individual subregions as shown 
in Figures [s] and 13-16, so that the fits are much better constrained than may be apparent from the combined distributions here. 
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Fig. 7. — Radial velocity vs. projected radius of the most 
negative- velocity stars in the NEl, NE2, and NE3 regions (red, 
blue, and magenta, respectively). Overplotted are data from two 
fields centered on the GSS south of M31 (turquoise triangles) 
( [Gilbert et al.| ||2009). The clusters of these turquoise triangles 
around v ~ —500 km s"-*^ and v ~ —390 km s"-*^ are the GSS 
and the secondary stream, respectively. The black crosses show 
six p oints for the predicted velocity of the NE Shelf (Fardal et alT] 
|2012| . The GSS appears to be present in the NE fields as slight 
concentrations around the black crosses. 



36 km for the thick and thin stellar disks, respec- 
tively. We expect our dispersion to be larger than the 
local value, primarily because the finite spatial extent of 
our subregions necessarily smears out the velocity dis- 
tribution. For example, if the spread in mean veloc- 
ity due to APA is ~ 20 km s~^ in a subregion, and 
the true velocity dispersion is ~ 40 km s~^, then 
we would expect to measure an effective dispersion of 
\/202 + 40^ = 45 km s~^. This effect is accentuated if 
there is an additional spread in mean velocity with ra- 
dius. Therefore, the fact that our cold component looks 
like the thick disk of [Collins et a l.|(2Qll) does not imply 
that the thick component dominates the stellar disk. 

A possible concern is that contributions from multiple 
stellar disk components (a combination of thin, thick, 
or extended) may invalidate our assumption of a locally 
cold disk. If this were the case, we would expect to see 
non-Gaussianity in the velocity signature of the cold com- 
ponent. However, Table [4] shows that that the velocity 
distributions of the cold components in most of the sub- 
regions are well fit by a single Gaussian. While this ob- 
servation does not say anything about the kinematical 
structure of the disk, it does suggest that our simple as- 
sumption is adequate for our purpose of describing the 
spheroid. 

It is also possible that a kinematically warm thick disk 
component would be incorporated into our Gaussian rep- 
resentation of the spheroid; however, the high goodness- 
of-fit statistic again suggests that this is not a signif- 
icant effect. In any case, the bias on our kinematical 
spheroid parameters induced by thick disk contamination 
confirms, rather than invalidates, our qualitative results, 
as discussed in § |4.5[ 

A closer look at the best-fit cold and hot components 
by subregion in Figures[8]and 13-16 show that the trends 
in the cold component with position angle match those 



expected for an inclined rotating disk. For example, in 
Figure § the mean velocity of the cold component (blue 
and green curves) transitions from —100 km s~^ in sub- 
region NEli, along the major axis, to —200 km s~^ in 
subregion NEl 7, farthest from the major axis. In other 
words, the absolute value of the offset of the mean ve- 
locity of the cold component from the systemic velocity 
of M31 moves closer to zero as we march away from the 
major axis. This progression is clear in Figures [13] and [M] 
for regions NE2 and NE3 as well, although it is less pro- 
nounced because these regions subtend a smaller range 
in PA. 

4.3. Effect of Tidal Debris Associated with the GSS on 
Spheroid Kinematics 

Accounting for the GSS and its associated debris has a 
significant effect on the measured kinematical parameters 
of the underlying s mooth inn er spheroid, supporting the 
observation of Fardal et al.l ([2007) th at the debris was 
visible in th e san iples oT TEata et ajj^ (|20C 



et al. 



(2006), and 



Merrett 



et al. 



Chapman 



I t is unclear whether the secondary str eam identified 



by K alirai et al. ( 2006b ) and confirmed by iGilbert e t al.| 
(|2Q09p is present in the NE1-NE3 regions. Figure [7[re^ 
veals a second clump of objects in the NE2 region, onset 
from the primary stream by 75-100 km s~^, the same 
separation as that between the two streams south of the 
galaxy. However, such clumps are barely, if at all, visible 
in the NEl and NE3 regions, and their exclusion does 
not significantly affect the kinematical parameters of the 
spheroid. Further observations are necessary to deter- 
mine if there is a second substructure and, if so, whether 
it has a physical connection to the GSS. 

Is the GSS the only source of nonvirialized substruc- 
ture biasing our measurements of the mean velocity and 
dispersion of the inner spheroid? While we cannot prove 
that all stars except those in the excised velocity ranges 
belong to either a disk or a perfectly smooth spheroid, we 
can show that the effect of other nonvirialized substruc- 
ture on our kinematical characterization of the underly- 
ing smooth spheriod is negligible. In star- co unt maps of 
the in ner regions of M31, such as those in |Ibata et al.| 
(2001), the GSS is by far the most prominent substruc- 
ture. Even so, GSS stars account for only a small fraction 
of our spectroscopic sample (see pink shaded regions in 
Figure [6]), and have a relatively small effect on our re- 
sults. Other, less-prominent substructures would (1) be 
nearly impossible to detect and account for, and (2) have 
a negligible effect on the measured kinematical parame- 
ters of the inner spheroid. 

4.4. Spheroid/Disk Membership Probability and 
Extreme Velocity Stars 

We can apply our subregion fits to quantify the 
spheroid membership probability for any star. At the 
location and velocity of each star in our sample, we calcu- 
late the ratio of the values of the best-fit hot component 
in that region and the disk and M32 components in that 
subregion. The disk and M32 membership likelihoods are 
calculated in a similar fashion. The results are sorted into 
three categories as follows: stars that are at least three 
times as likel y to be disk members as anything else (yel- 
low in Figure [11]) ; stars at least three times as likely to be 
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Fig. 8. — MCMC fits of kinematically hot (red) and cold (blue to green) components to each of 7 subregions in the NEl (inner northeast 
major axis) region plotted over velocity histograms of stars in each subregion. Subregions progress from NEli, straddling the NE major 
axis, to NEI7, farthest from the major axis. One curve is drawn for each component for each of the 32 walkers in the MCMC. Each 
violet distribution is the sum of the hot and cold distributions corresponding to a single walker. Velocity ranges excluded due to possible 
contamination by tidal debris from the GSS are shown in two ways: by the gray shaded regions, and by the stars shaded red on the 
histogram. The bottom middle panel shows the best-fit hot component (red) that, when added to the cumulative region cold distribution 
(blue), best fits the observed velocity distribution. The bottom right panel shows the final position of the walkers in parameter space. 
Ellipses show the mean and uncertainty of each of the parameters (v^ay) for each of the kinematical subcomponents. The dotted line in 
each panel marks the systemic velocity of M31 relative to the MW. 
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spheroid members as anything else (magenta); and other 
objects, including likely M32 members (green). As ex- 
pected, objects on the major axis are much more likely 
to be disk members than are those on the minor axis. 
Most important, we see likely spheroid members at all 
radii covere d by our sample . 



spheroidal systems. The bulges of spiral galaxies gener- 
ally fa ll on or above the line for rotationally-flattened sys - 



Recently, Caldwell et al. (2010) reported discovery of 
an "extreme velocity" star at a projected radius of 4 kpc 
(0.3°) along the SW major axis of M31. The star has 
a velocity of —780 km s~^, essentially excluding it from 
membership in the thin, cold stellar disk. Those authors 
attribute the star's highly negative velocity to possible 
membership in the GSS, even though it would have to 
be a 6(7 outlier of the stream velocity distribution. Our 
probability map demonstrates that this object can more 
easily be interpreted as a member of the spheroid, even 
at these large radii: it may be a 4cr outlier of the spheroid 
distribution. 

It is tempting to interpret this distribution of proba- 
bilities as a map of bulge-to-disk fraction. However, this 
statistic can be more reliably constrained using photo- 
metric light-profile fitting in conjunction with kinemati- 
cal decomposition. This study will be presented in Dor- 
man et al. (2012, in prep). 

4.5. Anisotropy 

We investigate the degree to which the flattening of 
the spheroid may be due to rotation by comparing the 
spheroid ellipticity e = 1 — b/a to the ratio v/ay. The 
value of e is uncertain at these radii. IPritchet fc vani 
|den Bergh (1994) mea sure 6 = . 4 with limited data at 
^proj = 10 kpc. Courtea u et al.| (|2011) perform a fit to 
a more extended data set to obtam values between 0.21 
and 0.37, depending on their bulge/disk decomposition 
and modeling technique, for a relatively small bulge with 
scale length '^1.0 kpc. 

Despite this range of possible ellipticities, v/ay of the 
spheroid in every region is lower than that of a rota- 
tionally flattened oblate isotropic rotator (black line in 
Figure 12). We measure v/ay = 0.23-0.52, but v^ot/cry 



of 0.87 or 0.54 would be required for rotation to flatten 
the spheroid to an ellipticity of 0.4 or 0.21, respectively. 
Hence, rotation alone probably does not account for all 
of the flattening of the spheroid. 

It is possible that anisotropy in the velocity ellipsoid 
can provide the remainder of the flattening. Anisotropy 
can be parameterized by 



1 



H. 



(6) 



where z points along the axis of symmetry of the 
spheroid, x is any orthogonal direction (say, r), and Ukk 
represents the pressure from velocity dispersion along di- 
rection k (e.g. Cappellari et al.||2007 ). In a system with 
an axisymmetric velocity ellipsoid, d = corresponds to 
a system whose flattening is unaffected by anisotropy, 
while S ^ 1 corresponds to a system whose flattening 
is almost entirely d ue to anisotropy. Using the iterative 
method described in Cappellari et ah] ([2007), we find that 
an anisotropy of S = 0.05-0.27 is required to explain a 
spheroid ellipticity of 0.21-0.4, given the mean value of 
v/ay in regions NE1-NE3 and SSW. 
We compare the inner spheroid of M31 to other 



tems fKor mendy fc iningworth||1982l |KQrmendy||1982b[ ) 

In contrast, so-called ''last-rotating ellipticals tend to 
lie between this line and the line S ^ 0.7eintr (magenta 



in Figure 12) ( j Cappellari et ah 2007). However, these 
comparison systems are observed at around one effec- 
tive radius (reff), whereas our kinematical measurements 
range from 1.5-14 reff. The inner spheroid of M31, then, 
more closely resembles the inner reff of a fast-rotating 
elliptical than the central bulge of a spiral galaxy. 



As suggested earlier in §4.2, it is possible that the 



spheroid velocity distribution may be contaminated by 
thick disk stars. If so, the true spheroid mean velocity is 
lower than we report, and the true spheroid dispersion 
higher. However, note that the v/ay of the spheroid is al- 
ready small enough to look more like an elliptical galaxy 
than the bulge of a spiral galaxy. The possible bias in- 
duced by a thick disk component would merely increase 
this effect, confirming our conclusion that M31's inner 
spheroid rotates unusually slowly. 

4.6. Effect of Target Selection Criteria on Velocity 
Distribution 

The majority of our spectroscopic targets were chosen 
on the basis of i^ magnitude only, but a small fraction of 
the targets in the MCT region were identified based on 
position in the PHAT CMD. This latter category con- 
tains proportionally more of the rarer metal-poor and 
metal-intermediate RGB populations. If metal-poor and 
metal-intermediate RGBs preferentially trace the kine- 
matically hot population, then inclusion of the PHAT- 
selected objects could bias our measurements, especially 
our spheroid membership probabilities. To test for bias, 
we performed a Komolgorov-Smirnov (K-S) test compar- 
ing the velocity distributions of the PHAT-selected tar- 
gets and the magnitude-selected targets in each subre- 
gion. The test confirmed that the two distributions were 
indistinguishable. 

Similarly, we used a K-S test to confirm that inclu- 
sion of crowded objects does not bias us towards one 
structural subcomponent. We created two lists of ob- 
jects: "crowded" (those sharing a slit with at least one 
serendipitously detected neighbor) and "isolated" (those 
without any such neighbors). Again, there was no singif- 
icant difference between the velocity distributions of the 
two categories in any subregion. 

5. SUMMARY &; CONCLUSIONS 

We have measured reliable radial velocities of over five 
thousand red giant branch stars in the inner 20 kpc of 
M31 using the Keck/DEIMOS multiobject spectrograph. 
Targets were selected using a series of statistical tests 
designed to identify isolated M31 members bright enough 
to yield quality spectra. 

By fitting a locally cold disk and a kinematically hot 
spheroid to the velocity distribution with an MCMC al- 
gorithm, we have measured the most probable kinemat- 
ical parameters v and ay of the red giant stellar popula- 
tion of the inner spheroid of M31 in each of five spatial 
bins. We find that, though the raw values are incon- 
sistent with a physical rotation pattern, accounting for 
the presence of tidal debris due to the GSS allows us to 
detect a significant spheroid rotation velocity. 
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TABLE 2 

Uncorrected Kinematical Spheroid Parameters 



Region 




in 


^proj (kpc) 


tteff (kpc) 


-^region 


V-VM31 (km s-i) 


av (km s ^) 


-^spheroid /-^total 


prob 


NEl 


0.40 


0.31 


6.45 


6.81 


1615 


+39.3 ± 13.9 


136.7 ± 5.4 


0.288 


0.731 


NE2 


0.48 


0.55 


10.11 


10.31 


1180 


+2.8 ± 24.6 


150.1 ± 12.0 


0.120 


0.136 


NE3 


0.64 


0.79 


13.93 


14.00 


859 


+48.0 ± 24.0 


131.4 ± 9.9 


0.203 


0.930 


SE 


0.26 


-0.14 


4.30 


6.91 


684 


+8.7 ± 10.3 


174.1 ± 10.5 


0.522 


0.121 


ssw 


-0.09 


-0.40 


5.69 


6.73 


1313 


-57.3 ± 8.9 


145.5 ± 5.5 


0.465 


1.5x10-^ 



TABLE 3 

Kinematical Spheroid Parameters Corrected for Tidal Debris 



Region 


vn 


^ (°) 


^proj (kpc) 


tteff (kpc) 


N 

-'^region 


V-VM31 (km s-i) 


ay (km s ^) 


-^spheroid /-^total 


prob 


NEl 


0.40 


0.31 


6.45 


6.81 


1615 


+41.8 ± 13.3 


134.4 ± 5.5 


0.297 


0.945 


NE2 


0.48 


0.55 


10.11 


10.31 


1180 


+31.3 ± 24.5 


135.3 ± 12.5 


0.130 


1.000 


NE3 


0.64 


0.79 


13.93 


14.00 


859 


+61.2 ± 16.44 


117.5 ± 8.8 


0.198 


0.998 


SE 


0.26 


-0.14 


4.30 


6.91 


684 


+9.0 ± 12.2 


144.5 ± 5.5 


0.571 


0.819 


SSW 


-0.09 


-0.40 


5.69 


6.73 


1313 


-57.1 ± 8.2 


145.4 ± 6.5 


0.484 


2.1x10-"^ 



We find that probable spheroid members are present at 
ah radii in our sample. When used in conjunction with 
integrated- light measurements, these membership prob- 
abilities will be a powerful tool for rigorous bulge/disk 
decomposition in the inner parts of M31. 

We also compare the v/(jy of the inner spheroid to 
those of other spheroidal structures. We find that rota- 
tion is insufficient to explain the ffattening of the inner 
spheroid; it more closely resembles an anisotropic fast- 
rotating elliptical than the bulge of a spiral or barred 
spiral galaxy as measured at ~ 1 reff. 

Our magnitude-limited survey biases us towards the 
bright, old stellar population. As more PHAT data be- 
comes available over the next few years, we will be able 
to select spectroscopic targets that represent a range of 
stellar populations within our magnitude range. We plan 
to present an analysis of the kinematics of different stel- 
lar populations in an upcoming paper. 
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Fig. 9. — Velocity (left) and dispersion (right) profiles for M31's inner spheroid without accounting for the GSS or associated tidal debris. 
Left: Mean velocity is plotted against the projection along the major axis coordinate. The SE minor axis (brown open square) is close to 
the minor axis, so we do not expect to measure a significant \v — Vsys\ due to rotation. The black line in the left-hand plot is the mean 
value of the four off-minor-axis velocity measurements. Right: Dispersion is plotted versus o-eff , the effective major axis coordinate for a 
spheroid with a 5:3 axis ratio. The best fit line to our dispersion measurements is shown by the solid black line, with the region within 
zb la of the best-fit line shaded in gray. Also plotted are limits on the bulge mean rotation velocity and velocity dispersion from Saglia] 
^ al. (201 0J black p oints with arrows), and velocity dispersion measurements further out in the spheroid by Chapma n et al.| (|2006) and 
[Oiibert et ^ ( |2"007| horizontal purple line). The orange dotted line represents a model dispersion profile that mcludes isotropic Dulge and 
Jialo stellar components, as well as a stellar disk and dark matter halo ( [Fardal et al.|2012| ). 




Fig. 10. 



Same as Figure |9] excluding stars in the velocity ranges corresponding to the GSS and associated tidal debris. 
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Fig. 11. — Probable subcomponent membership of each star 
based on location and velocity. Stars at least 3 times more likely 
to be a disk member than anything else are shown in yellow; those 
at least 3 times more likely than anything else to be a spheroid 
member are shown in pink; all other objects, including likely M32 
members, are shown in green. Region divisions are shown as solid 
lines. 
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Fig. 12. — Anisotropy diagram for spheroidal systems. The 
large filled squares show the ratio of the observed mean veloc- 
ity and the velocity dispersion against the intrinsic ellipticity of 
M31's spheroid. Error bars represent the range of reported el- 
lipticities converted to edge-on intrinsic ellipticities as explained 
in Cappellari et al. (2007). The black curve shows the expected 
value of the ratio v/ctv as a function of ellipticity for an isotropic, 
rotating galaxy. The magenta line approximates the behavior of 
an oblate anisotropic galaxy with anisotropy parameter 5 ~ 0.7e 
([Cappellari et al. 2007). Also shown are fast-rotating E/SO galaxies 
tem ISA U RUIN (turquoise squares) and bulges of barred [violet cir- 
cles, (Kormendy 1982b)] and unbarred [orange stars, ( Korm endy] 
&; Illingworth 1982)] spirals. The inner spheroid of M31 appears 
to be tlattened by velocity anisotropy in addition to rotation, and 
therefore more closely resembles an elliptical galaxy than a spiral 
galaxy bulge. 
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APPENDIX 

MCMC RESULTS FOR REGIONS NE2, NE3, SE, AND SSW 

Figures 13-16 are the analogs of Figure |8] for regions NE2, NE3, SE, and SSW, respectively. In this appendix, we 
briefly describe the results in each of these regions. 

NE2 Region 

Eleven stars are excluded from the MCMC fit in the NE2 region due to possible tidal debris membership; as seen in 
Figure [Tsj all are from the four subregions closest to the major axis. Exclusion drastically improves the chi-squared 
probability that the spheroid plus disk Gaussians are a good representation of the velocity distribution, raising it 
from 0.136 to 1.000. Exclusion also increases the mean velocity and decreases the velocity dispersion of the spheroid, 
bringing them much more into agreement with those measured in the other regions. 

Several trends characterize the subregion panels. First, the fraction of stars in the spheroid (parameterized by ratio 
of the areas of the red spheroid and blue/green disk curves) increases with increasing distance from the major axis. 
Second, the mean velocity of the cold component decreases with increasing distance from the major axis. Reading 
off the velocity axis of the parameter-space view in the bottom right panel, we see that the mean velocity of the cold 
component moves from about —90 to —170 km s~^. The minimum velocity is less negative than that in region NEl 
because region NE2 subtends a smaller range in PA. 

NE3 Region 

Nine stars are excluded from region NE3 as possible tidal debris contaminants. Exclusion slightly improves the 
reduced chi-squared of the fit. It also slightly increases the mean velocity of the hot component and decreases the 
velocity dispersion, bringing them more into agreement with those measured in the other regions. 

Because this region is centered farthest from the center of M31, it subtends the smallest range in PA of the three NE 
regions, and so the variation in cold component mean velocity between subregions is small. However, the parameter- 



space view in Figure 14 shows that the cold component kinematical parameters are very well defined. 



SE Region 

We exclude 36 stars from the tails of the velocity distribution in the SE minor axis region. While this does not 
skew the results of the MCMC fits in one direction or another, it does serve to increase the uncertainty in the final 
kinematical parameters; this effect can be seen in the large scatter in the walker positions in the parameter- space view 
in Figure [15] 

The minor axis is a saddle point in the velocity field of M31. Here, the observed rotation velocity changes from 
less negative than —300 km s~^ on the north side to more negative than —300 km s~^ on the south side. The mean 
velocities of the cold component reflect this transition. The parameter-space view in the final panel of Figure [15] shows 
that in subregion SEi (blue), which lies parallel to and just north of the minor axis, we measure a mean velocity 
slightly less negative than the systemic velocity of M31. In subregions SE2 (turquoise) and SE^ (green), situated 
south and north of the minor axis, respectively, we measure mean velocities more negative and less negative than 
—300 km s~^. Finally, we notice that the fraction of spheroid members is smallest in subregion SEi^ which is centered 
farthest from the center of the galaxy than the other two subregions. 

SSW Region 

Eleven stars fall into the velocity range excluded due to stream contamination in the SSW region. The exclusion 
does not have a significant effect on the mean velocity or dispersion of the hot component. The chi-squared probability 
of a good fit is lower in the SSW region than in the others, but exclusion of the possible stream stars improves the 
probability by a factor of about 15. 

Subregion SSW4 completely contains the galaxy M32, which is treated as a second cold component. (In a previous 
trial run, we allowed for the possibility of a second cold component in each of the five subregions, but the best-fit 
fraction of M32 stars in subregions SSWi_3^5 was zero.) The parameter-space view in Figure 16 shows that the mean 
velocity and dispersion of M32 are very well constrained around —200 and 27 km s~^, respectively. 

As in the NE regions, the fraction of stars in the spheroid is higher in the subregions farther from the major axis. 




Fig. 13. 



— Same as Figure [s] but for the NE2 region. 





V (km/s) V (km/s) 

Fig. 15. — Same as Figure[8] but for the SE region. SEi, SE2 and SE3 are the outer, inner south, and inner north subregions, respectively. 
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Fig. 16. Same as Figure [s] but for the SSW region. 
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Fig. 17. — Kinematical parameters of the cold component in each subregion, with la error bars, before exclusion of velocity ranges 
corresponding to GSS debris (black points) and after exclusion (colored points). Note that exclusion of tidal debris has a negligible effect 
on the cold component. These cold components correspond to M31's stellar disk. The left panel shows the rotation velocity of each cold 
component (computed using Equation|3]) vs. the mean rdeproj of each subregion. The three SE subregions have been excluded because they 
do not yield meaningful constraints on the rotation velocity. The right panel shows the velocity dispersion. The velocity and dispersion 
follow the trends expected for M31's stellar disk, supporting the idea that that our minimal assumption of a locally cold disk adequately 
accounts for disk contamination. 
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TABLE 4 

KiNEMATICAL DiSK PARAMETERS CORRECTED FOR TiDAL DeBRIS 



Subregion 




vn 


^proj (kpc) 


^disk (kpc) 


-^region 


V-VM31 (km s-i) 


av (km s ^) 


^cold/^total 


prob 


NEli 


0.284 


0.363 


6.311 


6.35 


220 


+192.4 ± 5.8 


56.4 ± 5.4 


0.725 


0.706 


NEI2 


0.287 


0.364 


6.36 


6.61 


227 


+186.8 ± 6.1 


64.2 ± 4.4 


0.806 


0.753 


NEI3 


0.286 


0.359 


6.33 


7.05 


236 


+179.4 ± 6.0 


58.6 ± 5.3 


0.729 


0.191 


NEI4 


0.321 


0.349 


6.56 


8.06 


234 


+184.5 ± 8.9 


63.8 ± 6.7 


0.602 


0.220 


NEI5 


0.362 


0.328 


6.74 


9.22 


164 


+181.8 ± 6.4 


47.0 ± 5.5 


0.599 


0.933 


NEle 


0.394 


0.281 


6.67 


11.01 


252 


+137.0 ± 6.0 


60.2 ± 6.5 


0.700 


0.862 


NEI7 


0.422 


0.180 


6.30 


14.64 


282 


+103.6 ± 5.2 


60.0 ± 6.8 


0.733 


0.520 


NE2i 


0.463 


0.597 


10.36 


10.42 


306 


+213.9 ± 3.7 


56.1 ± 2.8 


0.892 


0.767 


NE22 


0.459 


0.594 


10.30 


10.71 


264 


+206.9 ± 4.2 


52.3 ± 3.3 


0.943 


1.000 


NE23 


0.449 


0.583 


10.15 


11.30 


247 


+188.4 ± 4.7 


62.5 ± 4.0 


0.893 


0.005 


NE24 


0.481 


0.562 


10.25 


12.48 


135 


+198.5 ± 6.3 


53.5 ± 5.6 


0.749 


0.979 


NE25 


0.551 


0.434 


9.63 


14.17 


108 


+179.0 ± 8.7 


58.3 ± 6.6 


0.743 


0.987 


NE26 


0.568 


0.366 


9.26 


16.17 


120 


+133.5 ± 12.1 


80.1 ± 9.2 


0.861 


1.000 


NE3i 


0.653 


0.837 


14.56 


14.65 


339 


+206.6 ± 8.2 


47.7 ± 6.0 


0.770 


1.000 


NE32 


0.639 


0.805 


14.12 


14.68 


280 


+218.7 ± 9.4 


51.1 ± 7.1 


0.827 


0.104 


NE33 


0.609 


0.704 


12.83 


14.51 


240 


+206.2 ± 7.9 


54.6 ± 9.2 


0.820 


0.999 


SEi 


0.426 


-0.241 


6.68 


29.38 


130 


+16.0 ± 14.0 


52.8 ± 26.1 


0.448 


0.424 


SE2 


0.193 


-0.170 


3.62 


15.61 


355 


-25.5 ± 8.9 


49.8 ± 13.2 


0.345 


0.340 


SEs 


0.281 


-0.036 


3.95 


15.00 


199 


+46.3 ± 14.4 


77.9 ± 16.2 


0.567 


0.912 


SSWi 


-0.262 


-0.358 


6.07 


6.38 


254 


-200.9 ± 10.1 


74.2 ± 8.2 


0.647 


0.489 


SSW2 


-0.219 


-0.425 


6.54 


8.34 


149 


-183.7 ± 9.8 


74.5 ± 8.0 


0.779 


0.998 


SSW3 


-0.111 


-0.371 


5.31 


9.84 


175 


-145.3 ± 12.5 


61.6 ± 10.8 


0.449 


0.433 


SSW4 


-0.011 


-0.404 


5.54 


15.18 


608 


-88.5 ± 5.7 


33.6 ± 6.4 


0.191 


0.002 


SSW5 


-0.098 


-0.366 


5.19 


18.58 


127 


-54.7 ± 21.6 


61.8 ± 21.5 


0.340 


5.34 xlO-9 


SSW4 (M32)^ 


-0.011 


-0.404 


5.54 


15.18 


608 


+102.2 ± 2.7 


26.4 ± 3.2 


0.270 


0.002 



^ This second cold component in the SSW4 subregion corresponds to members of the M31 satellite M32. It is unrelated to M31's stellar disk and is 
only included in this table for the sake of completeness. 



